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We present a method to estimate the matrix element of the singlet axial current within 
a polarized proton state using lattice QCD. The method relies on using the Adler-Bell- 
Jackiw anomaly and gives the desired result in the chiral limit. We show that this method 
fails in the quenched approximation. For heavy quarks one does not expect much difference 
between simulations including dynamical quarks and those done in the quenched approx- 
imation. For that reason we explore numerical methods on an existing set of quenched 
configurations. The data obtained in this exploratory study show a poor statistical signal. 
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1. Introduction 

Four years ago, the European Muon Collaboration (EMC) measured the spin structure 
of protons using deep inelastic muon scattering from protons. Their principal result, which 
stirred a great deal of interest, was [1] 



»i 



/ gf(x)dx = 0.126 ±0.010 ±0.015 (1.1) 
Jo 



The quoted systematic error of 0.015 includes an estimate of the uncertainty due to the 
need to extrapolate the measured value of the structure functions to the regions above 
Bjorken x ~ 0.7 and below x ~ 0.02. The first moment of g\ is related to the spin 
structure of the proton through 



2 t g[(x, Q 2 )dx = (1 - ^±^)t[AAu(Q 2 ) + Ad(Q 2 ) + As(Q 2 )} + ... (1.2) 
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where corrections due to heavy quarks and higher twist effects have been ignored. The 
quantity Aq is the forward matrix- element of the quark axial current inside a proton, 



s^Aq{n 2 ) = (p,s\qi^^ 5 q 



\p,s). (1.3) 



Here is the covariant spin vector of the proton and a 2 is the renormalization point for 
the axial current. Since the singlet axial current picks up a small anomalous dimension 
only at two loops [2] we, therefore, ignore this dependence on the renormalization point 
in the following discussion. The EMC result (l.l)gives one linear combination of Au, Ad 
and As. A second relation is obtained from the neutron /3-decay 

(Au - Ad) = g A (1.4) 

where g A = 1.2573 ± 0.0028 is the axial vector coupling [3]. A third relation between the 
quark spin fractions can only be derived assuming SU(3) flavor symmetry. It is 

(Au + Ad- 2As) = 3F A - D A (1.5) 

where F A and D A are reduced matrix elements of the axial current. From an analysis of 
the strangeness changing semi-leptonic decays of hyperons, Jaffe and Manohar obtain [4] 

F A = 0.47 ± 0.04; ^ = 0.8^0.03. (1.6) 
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Combining Eqs. (1.1) - (1.6) they concluded that 

Au = +0.74 ±0.10 

Ad = -0.54 ±0.10 (1.7) 
As = -0.20 ±0.11 

and that the total axial charge of the nucleon is 

AS = Au ± Ad + As = 0.01 ± 0.29 (1.8) 

at Q 2 = 

Since in the parton model AS gives the fraction of spin carried by the quarks, the 
import of Eq. (1.8) is that the quarks contribute only a small fraction of the proton's spin. 
While this interpretation cannot strictly be derived from QCD, both this and the second 
consequence of Eq. (1.7) that there is a significant contribution from the strange quarks 
in the sea, are certainly challenges to our understanding of the structure of hadrons. 

The experimental situation has changed considerably in the last year. The SLAC 
E-142 collaboration has extracted AS and As from polarized electron scattering from a 
polarized 3 He target [5]. Their results are 

AS = ±0.57 ±0.11 

(1.9) 

As = -0.01 ±0.06 . 
More recently, the SMC collaboration [6] have reported the values 

AS = ±0.36 ±0.21 

(1.10) 

As = -0.07 ±0.06 

from polarized muon scattering off polarized deuterium. A global analysis of the data has 
also been made by Ellis and Karliner [7] including higher twist effects to compensate for 
the different values of Q 2 at which the three experiments are carried out. Using the values 

Fa = 0.46 ± 0.01; D A = 0.79 ± 0.01 . (1.11) 

taken from Ref. [8] they report 

Au = ±0.80 ±0.04 

Ad = -0.46 ±0.04 (1.12) 
As = -0.13 ±0.04 
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and 



AS = Au + Ad + As = 0.22 ± 0.10 . 



(1.13) 



The difference between these estimates of AS reflects to some extent the systematic error 
in the analysis due to the uncertainty in the contribution of small x region to gi, the 
uncertainty in including higher twist effects, and the uncertainty in the values of Fa and 
Da- Given the large spread in the results it is clear that the ability to calculate the matrix 
element of the singlet axial current in a proton will improve our understanding of this issue 
considerably. 

Numerical simulations of lattice QCD, in principle, provide an exact non-perturbative 
approach to calculating the matrix elements of the singlet axial current inside the proton. 
In Section 2 we present a method to extract the singlet matrix element using the Adler- 
Bell-Jackiw anomaly. In Section 3 we describe the problems with using this method in 
the quenched approximation and show that the matrix element is expected to diverge in 
the chiral limit. We discuss the lattice implementation in Section 4 and the subleties of 
calculating the topological charge density in Section 5. Assuming that for quarks with 
mass m q ~ m s the effect of vacuum polarization is small, we attempt an expoloratory 
calculation using the quenched approximation to explore numerical methods for a possible 
signal. The lattice parameters used are given in section 6 and the analysis of the data, 
which show a poor statistical signal, is presented in section 7. Our final conclusions are 
contained in Section 8. 

2. Relation between the axial current and the U(l) anomaly 

There are two possible Wick contractions that contribute to the matrix elements of 
the flavor singlet axial current inside a proton. These are shown in Fig. 1 and labeled as 
the valence and disconnected graphs. Unfortunately, of these two the one measuring the 
correlation of the disconnected loop with the proton has a very poor signal. A prelimimary 
lattice estimate found statistical errors nearly twice the measured value of the strange 
quark spin fraction [9]. Since this amplitude is not a priori tiny, contributes for each quark 
flavor, and is responsible for all of (p, s\ 5*7^75 s \p, s), we exploit an alternate method which 
is especially suitable to extract the total quark axial charge of the nucleon [10]. In the 
chiral limit the anomaly condition relates the divergence of the flavor singlet axial current, 
A M , to the purely gluonic operator FF 



ATT 



F --e 

£ v — 2 . 




(2.1) 
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where Nf = 3 is the number of active flavors. Taking the matrix elements of these operators 
inside a proton state and transforming to momentum space gives the relation 

AS = (p, s\A^, s) SfM = N f -± lim -^4{p\ s\TrF^F^(q)\p, s) (2.2) 

Z7T q— >0 q ■ S 

where q = p — p' and the two independent directions of proton spin are taken to be 
s/\s\ = ±q/\q\. Thus, operationally to get the final result we have to calculate the matrix 
elements between states with different momenta, divide by the momentum transfer, and 
then extrapolate to zero momentum transfer. The strong coupling constant a s should be 
evaluated at the same Q 2 at which the experimental result has been extracted even though 
the anomalous dimension of a s /2irTrFF is zero. The reason is that a s /2irTrFF and 
mix at 1-loop [2], have the same Q 2 evolution and thus have to be defined at the same 
renormalization point. 

Most numerical calculations are at present done in the quenched approximation. 
Therefore, in the next section we analyze the problems that arise in that case. We find that 
due to the vanishing of the rf mass as m q — ¥ one cannot calculate the spin asymmetry 
from the matrix element of FF alone and the proposed method cannot be used to get the 
desired result. 

3. The quenched approximation 

One a priori expects that vacuum polarization effects play an important role in the 
matrix elements of the singlet axial current due to the different behavior of rj' in the 
quenched and full theories. In particular, since the rj' becomes massless in the chiral limit 
of the quenched theory, the result may be singular. We show that this is indeed true and 
that in the quenched approximation Eq. (2.2) cannot be used to extract an estimate for 
the desired matrix element. 

We start with the Fourier transform of the anomaly relation for the matrix element 
of the flavor singlet axial current 

i QfX {f, s\ A^q) | p, s) = 2m q (p',s\P\ p, s) + N f % {p s\ Tr FF \ p, s) (3.1) 

where P is the SU(3) singlet pseudoscalar density operator. The first two terms have 
both connected and disconnected diagrams while the FF term only contributes through a 
disconnected diagram. These terms have single and double poles at q 2 = m 2 , due to the 
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f]' propagator [11]. We illustrate these in Figs. 2a-e. It turns out that in the q 2 — > and 
m q limit the desired matrix element is finite even in the quenched theory; nevertheless 
it cannot be calculated from (p' , s\ Tr FF \p, s) alone. 

We parameterize each of the three matrix elements present in Eq. (3.1) in terms of 
form- factors as 

(p', s\ Apfa) \p,s) = uij^ 5 uGi - iq^u-fsuG^ , 

(p',s\P\p,s) =u l5 uG p , (3.2) 
(p',s\TrFF\p,s) =u-f 5 uG F . 
Then, with the spinor normalization 

u u - ^ + M ?) l + ^ (33) 

Up,sU p , s - 2 , [6.6) 

one has the relation G^(q 2 = 0) = AS. The singularities in these form factors with respect 
to q 2 and due to the rj' propagators are 



G F {q 2 )= (2 h 2 , +F. 

Note that the rf can couple to the proton through pseudoscalar and axial couplings. How- 
ever, both of these are proportional to q^u^u for on-shell matrix elements. The form factor 
Gf gets contributions from contact terms and the exchange of 1 + states, but not from rf 
propagating as an intermediate state. Thus, only G2 contains rf pole terms. Equating the 
single and double pole terms gives two relations 

m 2 ,a 2 = 2m q p 2 

m v> a i = — T^ + ^Pi+Nf^fi 
between the residues. Now using Eqs. (3.1) and (3.2) we can write 

2M P Gt{q 2 ) = -fGitf) + 2m q G p (q 2 ) + N f ^G F (q 2 ) (3.6) 
which has the limiting behavior 

lim 2M P Gf(q 2 ) = -a x + 2m q P + N f f-F, 

q 2 — ¥0 • 7r 

lim 2M P Gt{q 2 ) = -a, - q 2 G 2 + N f £F, 
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(3.7) 



showing that Gf is finite in the chiral and q 2 — ¥ limit. In the lattice calculation we first 
take q 2 — ¥ and then extrapolate to the chiral limit. In this double limit one gets the 
result 

2M P G?(q 2 = 0) = - ai +N f %F, 

2m 9 / p 2 \ , A7 a .,* /i x ( 3 - 8 ) 



The second expression is obtained using Eq. (3.5). Using the anomaly relation (2.2) to 
calculate the spin asymmetry gives only the term proportional to Nf^. This diverges 
in the chiral limit and there is no obvious way of extracting the physical answer from it. 
Thus, to get the result in the quenched approximation, the best option is to calculate the 
connected and disconnected diagrams for the axial current itself. 

In the full theory, there are no double poles as the diagrams in Figs. 2b-d, along with 
all diagrams related by the insertion of vacuum polarization bubbles, are part of the single 
r]' propagator shown in Fig. 2f. Therefore, keeping only single pole terms in Eq. (3.4) and 
repeating the above analysis, one finds the relation between the residues to be 

m^at =2m qPl +N f ^f 1 . (3.9) 

In the double limit q 2 — > 0, m q — > one gets 

2M P G£(q 2 = 0) = -ai + N f %F, 

h ^ ( 3 - 10 ) 

which is precisely the anomaly relation (2.2). 

To summarize, our analysis exposes a very serious shortcoming of the quenched ap- 
proximation. The matrix element one gets using Eq. (2.2) is expected to diverge in the 
chiral limit. The most one can hope for is that for quarks of mass > m 3 there is little effect 
of quark loops as one is far from the chiral limit. In this region one can use quenched lat- 
tices to develop numerical tricks to improve the signal that will be relevant for calculations 
in the full theory. Thus our present numerical work should be regarded as an exploration 
of methods. 
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4. Lattice implementation 

On the lattice we calculate the following correlation function 

Re (^j< N ir (t)(P+-P-)TrF flv F flv {T,$N p (0) )) (4.1) 

where Np(t) is a creation/annihilation operator for the proton at time t and momentum 
p, P± are the spin projection operators, and in Euclidian space P + — P_ = ^75 projects 
out the spin asymmetry. In order to extract the desired matrix element we have to remove 
the amplitude for creating and annihilating the proton and the exponential damping in 
time from the above correlator. In addition it is necessary to choose £>>r>>0so that 
only the lowest proton state contributes. To achieve this we calculate the following ratio 
of correlators 



Re ^Tr<%,(*)(P+ - P_)i^i^(r, g)%(0)) 

TZ(t, r) = —i , '- (4.2) 

Re(Tr<%(*)^%(0)}J 

It is easy to see that in this ratio the normalization of the proton field, including the 
Wuppertal smearing factors, cancel. In addition, the exponential decay factor e~ ENt is 
cancelled if the magnitude of the momentum is unchanged between the inital and final 
proton state, i.e. [p| = \p'\. 

We have calculated this ratio for three classes of momentum transfer: (1,0,0) — > 
(0,1,0), (1,0,0) ->- (-1,0,0), and (2,0,0) -> (-2,0,0). In each case the results are aver- 
aged over all permutations of the axes, and over forward and backward propagation on a 
given background configuration. Statistical errors are then computed by a single elimina- 
tion jackknife method regarding each of the 35 data points as statistically independent. 

The cases (1,0,0) -> (-1,0,0) and (2,0,0) ->• (-2,0,0) are simple to analyze. We 
choose the two independent directions of proton spin to be s/\s\ = ±p/\p\. The case 
(1,0,0) — > (0,1,0) has the advantage that, under the restriction that |p| = \p'\, the 
momentum transfer is the smallest possible. In this case we take the spin projection to be 
along the direction of momentum transfer q = (p — p 1 ). 



5. Topological density 

The calculation of the topological density Q(x) = TrF^F^ on the lattice has not 
proven to be straightforward. In the continuum limit, we require that the topological 
charge, given by Q = — 1/167T 2 J x Q(x), take on integer values and be related to fermionic 
zero modes, i.e. Q is the difference between the number of left and right handed zero 
modes. These properties hold on the lattice near the continuum limit where the lattice 
gauge fields are smooth. The problem with present lattice calculations (j3 « 6.0) is that 
the gauge fields are not smooth. As a result there exist lattice artifacts, called dislocations, 
that are caused by ultraviolet fluctuations, and these affect the calculation of Q and need 
to be removed before one can get sensible results. Also, on the lattice Q(x) mixes with 
the other dimension four pseudo-scalar operator mip^ip and only a linear combination is 
equal to the divergence of the axial current. Our method of using the Adler-Bell-Jackiw 
relation circumvents this problem, on the other hand the results are only valid in the chiral 
limit . 

We have used the naive lattice transcription of Tr F^F^ where the gluon field strength 
Ffxv{ x ) is defined as the traceless anti-hermitian part of the plaquette. To improve the 
estimate with respect to 0{a) corrections we use the symmetrized average of the four 
plaquettes in the [/i, u] plane [12]. While this field theoretic definition (which we shall call 
the plaquette method) has the correct continuum limit, it is well known that Q does not 
take on integer values [13] at finite lattice spacing and it is believed that the contribution 
of lattice artifacts may be large. We discuss this point and operator renormalization in 
some detail below. 

A number of numerical methods have been proposed to smooth out the gauge fields 
without changing the topological character of the configuration. These methods, which go 
under the name of "cooling" and "smearing" , have been compared in, for example, Ref. [14] 
but it is still not clear whether the smoothing process washes out only the lattice artifacts 
or in addition some physical features of Q(x). It is therefore important to calculate Q(x) a 
number of different ways and check whether the result for the matrix element is the same. 

The success of all the plaquette methods relies on the validity of a decomposition of 
the form [15] 

Q L (x) = (Z(/3) + C(x))Q(x)+ri(x) (5.1) 

where ((x) and r](x) are random variables with zero mean, and Z(j3) is the renormalization 
constant. Numerical data for the 0(3) non-linear sigma model in two dimensions show 
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that such a decomposition is justified and Z(j3) agrees with the 2-loop perturbative result 
[15]. We shall assume that this is true also for 4-D gauge theories. The problem then 
reduces to one of getting adequate statistics, i.e. averaging over a sufficient number of 
gauge configurations such that the correlations of ((x)Q(x) and r](x) with the proton 
correlator are zero. This assumption can be tested by comparing results using different 
methods of calculating Q(x). 

Assuming that the decomposition given in Eq. (5.1) works, one still needs the mul- 
tiplicative renormalization factor Z((3). This has been calculated in perturbation theory 
and the result of Ref. [16] is 

*</» = l + J + ! (5-2) 

where Z\ = —5.4508 and Z^ ~ 11.1. (In the 2-loop part Z^ only the dominant graphs have 
been evaluated). Thus Z(j3 = 6.0) ~ 0.4. The renormalization factor was also estimated 
using a cooling procedure [16], which gave Z((3 = 6.0) ~ 0.3 ± 0.07. 

An estimate of Z(f3) can also be obtained using the mean-field method of Lepage and 
Mackenzie [17] 

Z = Uq = Up laquette fa 0.35 (5.3) 

where Uq is the expectation value of the plaquette. The three estimates are consistent and 
show that the renormalization effects are very large. In this study we shall present results 
using the mean-field value given in Eq. (5.3). 

An alternate approach to the plaquette method is to use Luscher's geometrical def- 
inition of Q(x). This, by construction, guarantees that Q takes on integer values and 
has the correct continuum limit. However, it has been shown that even in this method 
dislocations contribute to Q(x) for the simple Wilson action [18]. The suggested cure is 
to use an improved gauge action which increases the action of dislocations so that they 
are supressed for [5 > 6.0 and vanish in the continuum limit. Recently, Schierholz has 
argued that at corresponding values of the lattice scale, dislocations are also suppressed in 
simulations with light dynamical fermions [19]. We are in the process of developing codes 
to calculate Q(x) using Luscher's method so as to compare the results with the plaquette 
method (with and without using smeared or cooled lattices). 

To summarize, the plaquette method for calculating Q(x) is computationally much 
faster, however it does assume that Eq. (5.1) is valid and that Z((3) can be estimated 
reliably. Since all lattice methods for calculating Q(x) have potential problems and the 
statistical signal is poor, we consider it worthwhile to compare results for the matrix 
element with different definitions of Q(x) on the same set of lattices. The use of the 
plaquette method should be considered as the first step towards this goal. 
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6. Details of the lattices and quark propagators 

The calculation uses 35 quenched gauge configurations of size 16 3 x 40 at = 6.0. The 
Wilson action quark propagators were calculated on doubled lattices (16 3 x 40 — > 16 3 x 80) 
using Wuppertal sources. The quark masses used are k = 0.154 and 0.155, corresponding 
to pions of mass 700 and 560 MeV respectively. More complete details of lattice generation 
and propagator inversion are given in Ref. [20]. These parameter values are far from the 
real world, and as discussed in Section 3 we do not expect realistic results in the quenched 
approximation. Nevertheless, quenched lattices form a good laboratory to test numerical 
methods, so the goal of our calculation is to explore techniques to enhance the signal. 
Unfortunately, our conclusion is that an improvement in statistics and numerical methods 
will be required to get a reliable signal in eventual simulations with dynamical quarks. 

7. Analysis 

There are three quantities which control the quality of the final signal. We discuss 
each of these below. 

7.1. The operator (FF): 

The expectation value of Q(t, q) should be zero on each time slice and for each value 
of momentum transfer. Our data satisfy this condition within the error estimate for all 
three cases of momentum transfer. The data for the best case, p = (1, 1,0), is shown in 
Fig. 3 and we find that the signal is very noisy and that the correlations between successive 
time- slices are large. 

7.2. Proton Correlators 

At large time separations the proton correlator should be saturated by the lowest 
state. The data for p = (1,0,0) are shown in Fig. 4. We find a plateau in the effective 
mass for p= (0,0,0)andp= (1, 0, 0) and a marginal signal for p = (2,0,0). The saturation 
(> 95%) by the lowest state sets in for t > 8 in each case as examplified in Fig. 4. This 
implies that in order to get a reliable estimate for the matrix element, the insertion of FF 
should be at least 8 time slices from both the proton source and sink. It is hoped that 
some of the effects of the higher states in the proton correlator cancel in the ratio defined 
in Eq. (4.2) so that one can extract reliable estimates of the ground state matrix element 
from slightly smaller separations. Unfortunately, as we discuss below the signal in our data 
is not good enough to draw any reliable conclusions about the stability of the results. 
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7.3. Correlation of FF with the Proton 

The correlation between FF and the proton occurs as a result of gluon exchange 
in the i-channel. This is shown schematically in Fig. 2e. The analogue of this in the s 
channel is the off-diagonal matrix element for the production of a pair of nucleons by a 
glueball. Based on glueball calculations one expects that to get a good signal will require a 
large number of configurations (many more than the 35 configurations used in this study). 
However the situation may not be this bad. The correlation of FF with the proton comes 
mainly from contact terms, i.e. correlated fluctuations of the same link variables that enter 
into FF and also into the Dirac equation for the quark propagators. This suggests that it 
is the short distance correlations which give the dominant contribution. So the signal one 
hopes for may have more in common with the computation of glueball propagators at zero 
and small separations than that at large separations as needed for extracting masses. It is 
this physical picture that gives us hope of a reasonable signal in the present calculation. 

The data for the three different cases of momentum transfer, (1,0,0) —¥ (—1,0,0), 
(2,0,0) — > (—2,0,0) and (1,0,0) — > (0,1,0), are shown in Tables 1 for the heavier mass 
k = 0.154. All numbers should be multiplied by 10~ 2 a s to get AS. The data are shown 
as a function of the time slice at which the proton is annihilated (Tp) and at which the 
operator (T^p) is inserted. The proton is created at time-slice T = 0. We find that the 
statistical errors in the results are such that no signal is evident, and that these data are 
statistically indistinguishable from zero. 

Qualitatively, we find as expected that the statistical errors are best for the case 
(1,0,0) — y (0,1,0), somewhat larger for (1,0,0) — > ( — 1,0,0) and still larger for the case 
(2, 0, 0) —¥ (—2, 0, 0). Also the errors increase with Tp and the maximum separation we can 
probe is Tp = 9, significantly short of the requirement that the insertion of FF should be 
at least 8 time-slices from either end. In view of this lack of signal, some future possibilities 
to improve the signal are to work on a larger lattice in order to reduce the value of \q\, to 
use alternate definitions of Q(x), to find better smearing for quark propagators so that the 
lowest state saturates the correlation function at short times, and to significantly increase 
the statistics. 
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8. Conclusions 

We have presented a method for measuring the spin asymmetry of the proton using 
the Adler-Bell-Jackiw anomaly. Unfortunately, this method fails in the quenched approx- 
imation due to the r/' pole at q 2 = 0. We discuss the calculation of the topological charge 
density using the simplest field-theoretic method and explore numerical methods for ob- 
taining a signal for the matrix element of FF inside a proton using an existing set of 
quenched lattices. We do not find a reliable statistical signal and expect that future calcu- 
lations using unquenched lattices will be needed to explore alternate methods to improve 
the signal. 
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-3(11) 










_ 7 


28(22) 


-28(22) 


-23(19) 


-21(17) 


23(18) 


30(26) 






Tp 


= 7 


28(30) 


-20(41) 


56(46) 


-8(27) 


-66(36) 


-93(41) 






Tp 


= 7 


-2(16) 


2(17) 


19(18) 


22(14) 


-14(19) 


15(15) 






T 


Q 

= 


36(29) 


-29(29) 


-24(24) 


-17(29) 


28(27) 


49(37) 


40(36) 




T P 


= 8 


44(40) 


-6(59) 


63(67) 


1(50) 


-96(55) 


-207(75) 


-47(78) 




T P 


= 8 


-1(19) 


13(24) 


20(23) 


31(19) 


-13(25) 


26(23) 


-6(25) 




T P 


= 9 


52(35) 


-15(39) 


-25(33) 


-9(41) 


26(36) 


61(43) 


50(50) 


-22(40) 


T P 


= 9 


48(73) 


22(94) 


38(94) 


14(93) 


-134(95) 


-305(118) 


-95(102) 


-13(99) 


T P 


= 9 


-6(27) 


32(31) 


39(28) 


39(24) 


-21(29) 


34(32) 


-7(30) 


-21(25) 



Table 1: The data for the matrix element as a function of T-^p (time-slice of operator 

insertion) and Tp (sink time-slice of the proton) at k, = 0.154. The source time-slice of 

the proton is held fixed at T = 0. The three numbers at each value of Tp and T„ ^ are 

r r 

for the three cases of momentum transfer, i.e. q = (2,0,0), q = (4,0,0), and q = (1, 1,0) 
respectively. All numbers should be multiplied by 10 _2 a s to get AE. 
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A) 



B) 




Fig. 1: The connected (a) and disconnected diagram (b) that contribute to the matrix 
element of the singlet axial current in a proton. 
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Fig. 2: Diagrams showing the contribution of r)' poles. The symbol X denotes the insertion 
of the current and the dotted lines represent gluon exchange. Fig. 2a shows the single pole 
in the connected diagram, Figs. 2b-e show double and single poles in the disconnected 
diagram. In the full theory Figs. 2a-d are part of the same single pole shown in Fig. 2f. 
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Fig. 3: The expectation value of TrFF as a function of time slice for q = (1, 1, 0). Different 
symbols have been used to help distinguish between the positive and negative values. 
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Fig. 4: The proton correlator and the single mass fit. We also show the saturation by the 
lowest state and the effective mass plot as a function of the separation t. The data are for 
p= (1,0,0). 
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